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Abstract 

Densities from relativistic mean field calculations are applied to con- 
struct the optical potential and, hence calculate the endpoint of the rapid 
proton capture (rp) process. Mass values are taken from a new phe- 
nomenological mass formula. Endpoints are calculated for different temperature- 
density profiles of various X-ray bursters. We find that the rp process can 
produce significant quantities of nuclei upto around mass 95. Our results 
differ from existing works to some extent. 

1 Introduction 

Proton capture reactions at a very low temperature play an important role in 
nucleosynthesis process. Most importantly, in explosive nucleosynthesis (e.g. 
an X-ray burst), the rapid proton capture (rp) process is responsible for the 
production of proton-rich isotopes upto mass 100 region. In nature, the proton 
capture reactions, important for nucleosynthesis, usually involve certain nuclei 
as targets which are not available on earth or cannot be produced in terrestrial 
laboratories with our present day technology. Therefore, theory remains the 
sole guide to extract the physics. 

In our present work, we have studied the endpoint of the rp process in a 
microscopic approach using a new phenomenological mass formula 1 . In a 
similar work, Schatz et al. [2] calculated the endpoint of rp process using mass 
values from a Finite Range droplet Model (FRDM)[3] calculation and proton 
capture rates from Hauser-Feshbach code NON-SMOKER 4 , 5 . We will show 
that the results of our calculation are different from their observations. In the 
present work, we will concentrate only on X-ray burst scenarios, which have 
typical timescale of 100 seconds and a peak proton flux density of the order of 
10 6 gram I cm 3 . This type of burst provides a highly proton-rich environment 
around the peak temperatures 1-2 GK. We try to look at different models of 
the X-ray burster and find out the endpoint of the rp process nucleosynthesis. 



1 



2 Methodology 



When an X-ray burst takes place, a proton-rich high temperature environment, 
which triggers the rp process, is created. The process passes through nuclei near 
the proton drip line, not available on earth. In regions far from the stability 
valley, rates derived from phenomenological calculations may not represent the 
reality very well, leading to considerable uncertainty in the process. Very often, 
the reaction rates are varied by a large factor to study their effects. On the 
other hand, in a microscopic calculation, uncertainty in reaction rates can be 
reduced and therefore, this approach is expected to give a more accurate result 
for unknown mass regions. In a previous work[Bj, we have shown that the rates 
may vary at most by a factor less than two when the cross-section values range 
over four orders of magnitude. 

A microscopic calculation has been performed to evaluate proton capture 
rates for the nuclei involve in the rp process in the present work. We use the 
spherical optical model to calculate the rates of the relevant reactions. As most 
of the nuclei involved in the process lie around the drip line, experimental den- 
sity information are not available. Hence, theoretical density profiles have been 
calculated using relativistic mean field (RMF) theory. In the present work, we 
use the FSU Gold Lagrangian density [7'. and solve the RMF equations in the 
spherical approximation for the nuclei involved in the rp process. This La- 
grangian density, containing additional nonlinear terms for the vector isoscalar 
meson self interaction and an interaction between the isoscalar vector and the 
isovector vector mesons, has been found to be very useful in describing nuclear 
properties throughout the mass table [See e.g. Bhattacharya et a/[SJ H] and 
references therein]. 

The microscopic optical model potential for capture reactions are obtained 
using effective interactions derived from the nuclear matter calculation in local 
density approximation, i.e. by substituting the nuclear matter density with 
the density distribution of the finite nucleus. In the present work, we have 
constructed the potential by folding the density dependent M3Y fDDM3Y)|101 
ITT1 H2] interaction with densities from RMF approach. This interaction was ex- 
tracted from a finite-range energy- independent G- matrix element of the Reid po- 
tential by adding a zero-range energy- dependent pseudo-potential and a density- 
dependent factor. The interaction at the point f is, thus, given by 



The density dependent factor g(p) has been chosen of the form C(l — ftp 2 / 3 ) 
from the work by ChaudhurijTTJ [12] where the constants were obtained from a 
nuclear matter calculation as C = 2.07 and (3 = 1.624 fm 2 . We have used this 
form in our calculation without changing any of the above parameters. 




t M3Y = 7999 :_ _ 2134—— + J m (E)S(r) (2) 
4r 2.5r 



for r in fm, and Jqq(E) is the zero range pseudo potential, 




(3) 
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We have also fixed the various parameters and prescriptions in the Hauser- 
Feshbach calculation for the relevant mass region by comparing our results to 
the experimental low energy proton capture cross sections for these nuclei. Our 
method of calculation for mean field and proton capture rates has been described 
in our earlier works [6l [13] in detail. The computer code TALYS[T4] has been 
used for rate calculation. 

Binding energy of nuclear ground state is one of the most important inputs 
in the study of astrophysical reactions. Experimental measurements are very 
difficult to perform in nuclei far from the stability line. Therefore, one has to 
take recourse to theoretical predictions. Though we have used a mean field 
calculation to extract the nuclear density profiles, no RMF calculation has been 
able to achieve a prediction of mass values with accuracy sufficient for studying 
the proton drip line. In fact, even Skyrme Hartree-Fock calculations can predict 
the mass values with an root mean square (rms) error slightly less than 0.6 
MeV only. Thus, in the present work, we have obtained the mass values from a 
newly developed mass formula[l . It uses a purely phenomenological form with 
empirically fitted parameters and predicts the known mass values of 2140 nuclei 
with an rms error of 0.376 MeV. In view of the success of the formula to predict 
the proton dripline and rp process upto mass 80 region 15 and to predict the 
peaks in r process [TB] quite well, it will be interesting to see the effect of this 
mass formula on rp process beyond mass 80 region and to the endpoint of the 
rp process. 

In an X-ray burst environment, a nucleus (Z, A) may capture a proton to 
form the nucleus (Z + 1,A + 1) . However, this process has to compete with 
its inverse, i.e. photodisintcgration by emitting a proton at high temperature 
[(7,p) reaction]. A negative or a small positive value of the proton separation 
energy implies that the inverse reaction dominates and the rp process stalls at 
that point, the so-called waiting point. Therefore, only a two proton capture 
process can bridge the waiting point nuclei. The bridging mechanism has been 
discussed in standard text books (for example the book by IniadisfTT]). For the 
X(p,-f)Y reaction, the rate A of the inverse process (j,p) is related to the the 
proton capture rate by the reciprocity theorem and is of the form [17] 



in the unit of sec -1 for X(p,~/)Y process. The forward reaction rate, denoted 
by N (av)px^Y-y, is expressed in cm 3 mol~ 1 sec~~ 1 . The temperature T is in 
GK (10 9 K) and Q is the ground state Q-value of the (p, 7) reaction expressed 
in MeV. The normalized partition functions, Gx and Gy, have been obtained 
from Rauscher et al.{l& For protons, we use the standard values, G p — 1 and 
J p = 1/2. As evident from the above expression, the Q-value, appearing in the 
exponent, plays a vital role in the whole process. 

Apart from the above processes, a nucleus can decay by emitting beta par- 
ticles while, for higher mass isotopes, a-decay is another probable channel. In 
this work, the measured half life values for /3-decay have been taken from the 
compilation by Audi et o/.[Tl] except in the case of 65 As, which is taken from 
the experimental measurements by Lopez et aZ.|20j In absence of experimental 



A = 9.86851 x 10 9 T2 



1 ( M P M X \ 2 (2J P + 1) (2J X + 1) G p Gx 

2 V My ) (2J Y + 1) G Y 
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data, half life values have been taken from the calculation by Moller tt aL[5] 
both for j3- and a- decay. Taking into account all the above processes, we have 
constructed a network extended upto A=110 region to study the nucleosynthesis 
path and relative abundances of the elements at any instant. 

3 Results 

The rp process paths are shown in Fig. 1 and Fig. 2 for temperatures 1.2 and 
1.5 GK, respectively, using a constant proton flux density of 10 e gram/cm :i , a 
proton fraction of 0.7 and 100 seconds burst duration. Here, black lines indicate 
the path along which the major portion of the total flux flow whereas gray lines 
indicate the minor paths. The filled boxes in figures indicate the waiting points. 
As evident from these figures, rp process paths depend on the temperature of 
the environment to some extent. For example, at T= 1.2 GK, major portions 
of the rp process flux at the waiting point nucleus 64 Ge convert to 66 Se by two- 
proton capture process. Less than 1% of the total flux flows through the /3-decay 
channel and follow the paths showed by gray lines in Fig. 1. In contrast, at 
temperature T= 1.5 GK (in Fig. 2), the probability of two-proton capture of 
64 Ge gets reduced and the rp process path bifurcates from the waiting point 
nucleus almost in an equal proportion. This suggests that, at 1.2 GK, proton 
capture by 64 Ge dominates over its inverse process, i.e. photodisintegration. 
As the temperature increases, contributions of photodisintegration process in- 
creases and therefore a large fraction of the total abundance chooses another 
path through more stable nuclei such as 65 Ge, 66 As, etc. Similarly, near other 
waiting points (viz. 72 Kr, Sr etc.) the abundance flow pattern changes with 
changing temperature. It is also evident from above the figures that, beyond 
mass 80, temperature change can hardly affect the scenario. 

As shown in Fig. 1 and Fig. 2, above A=100 region, the rp process continues 
through proton capture by In isotopes and /3-decay of Sn isotopes. Here, 100 In 
captures a proton to form 101 Sn which completely decays to 101 In as 102 Sb is 
proton unbound. In turn, 101 In undergoes proton capture and further exhibits /3- 
decay. The process of proton capture followed by consecutive /3-decay continues 
and relative abundances of nuclei decrease as one proceeds towards higher mass 
region. Ultimately, less than 0.001% of the total flux can reach 106 In, according 
to our present calculation. 

According to the calculations of Schatz et al. 2J , a significant portion of the 
105 Sn captures a proton to form 106 Sb, as 106 Sb has a sufficient positive proton 
separation energy (0.59 MeV from FRDM[3] calculation). Another proton cap- 
ture leads to 107 Te, which instantly undergoes a-decay to 103 Sn. Thus, the rp 
process ends in the SnSbTe cycle. 

In contrast, our calculation does not go through the SnSbTe cycle. It is 
clear that the SnSbTe cycle may not occur under two scenarios. Firstly, if 106 Sb 
isotope be very loosely bound (proton separation energy 0.119 MeV according 
to the mass formula [I]), any 106 Sb that is formed by a proton capture instantly 
reverts back to 105 Sn. Second scenario occurs when the proton capture rates 
are too small for 105 Sn to initiate a proton capture process. 

In a previous work [15]. we have shown that a small fluctuation in the mass 
values of waiting point nuclei in mass 60-80 region may affect the effective half 
life and thus, proper knowledge of ground state binding energy is necessary to 
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understand the bridging phenomena of a waiting point nucleus below mass 80 
region. However, as we move towards the higher mass region, we find that small 
variations in binding energy do not affect the rp process path significantly. 
Taking the rms error into account for the proton separation energy of 106 Sb 
(0.119 ± 0.376)MeV, we have checked whether the proton capture on 105 Sn 
can dominate over its inverse process. We find that the fraction of the initial 
flux entering into the SnSbTe-cycle is negligibly small. We have repeated our 
calculation with the proton separation energy of 106 Sb isotope from a recent 
experiment by Elomma et al. [21 , viz. 0. 428(0. 008)MeV. We find that our results 
remain almost invariant. Repeating the entire calculation with ground state 
binding energies from FRDM[3 calculation or the Duflo-Zuker[22] mass formula 
do not alter this conclusion. Hence, we see that for a reasonable variation of 
mass values the rp process fails to enter into the SnSbTe-cycle for an X-ray 
burst of 100 seconds duration. If we consider an X-ray burst having duration 
greater than 150 seconds, a very small fraction of total flux, less than 0.001% 
of the original, enters the cycle. 

Schatz et ai.[2] have used capture rates different from ours. In Fig. 3, 
we have plotted proton capture rates with temperature for the isotopes 105 Sn 
and 106 Sb respectively. The different values have been indicated as follows: 
Pres - Present work and NON-SMOKER - NON-SMOKERgl [5] results used 
by Schatz et al. [2] In case of 105 Sn nuclei, it is evident from Fig. 3 that proton 
capture rates from our calculation are approximately 10 2 times smaller than 
the rates from NON-SMOKER calculation. It seems, the difference is mainly 
due to the fact that they have used the form of the interaction from Jeukcnnc 
et aL|23] which is different from our case. For better understanding, we have 
repeated the entire calculation using reaction rates from Hauser-Feshbach code 
NON-SMOKER[H [5]. In this case, considerable fraction of the initial flux (~ 
0.1%) enters into the SnSbTe cycle. Thus, we may conclude that the proton 
capture rates, and not the mass values, determine the endpoint of rp process. 
In our opinion, our approach in fixing the parameters in the reaction calculation 
by fitting known reaction rates, and extrapolating the calculation to unknown 
reactions, may be relied upon to predict the endpoint correctly. 

There are many X-ray burster models in literature with various density- 
temperature profiles. Our present goal is to study how different temperature- 
density profiles can affect the relative abundances of the nuclei produced by rp 
process and the endpoint. In Fig. 4, we have plotted the relative abundances 
of elements versus mass number at 100th second of the burst in various models, 
as described below. 

In the first case (Model-I) , we assume a constant density-temperature frame- 
work with temperature T = 1.5 GK, density p = 10 6 gram/cm 3 and pro- 
ton fraction value 0.7. Model-II[2 describes a situation where ignition takes 
place at a constant density of p = 1.1 x 10 6 gram/ cm 3 and the burst reaches 
a peak temperature T pea k =1.9 GK after 4 seconds while the cooling phase 
lasts for approximately 200 seconds. In another example, taken from the book 
by Illiadis[T7] (Model-Ill), nuclear burning starts with temperature and den- 
sity values of T = 0.4 GK and p — 10 6 gram/ cm 3 , respectively. After 4 sec- 
onds, a maximum temperature of Tpeak — 1-36 GK and a minimum density of 
Ppeak — 5 x Ingram/ cm 3 are achieved. After 100 seconds, the temperature 
drops to T = 0.7 GK and the density increases to p = 1.4 x 10 6 gram / cm 3 . For 
all the above models, it is assumed that the X-ray burst environment is suffi- 
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ciently proton-rich to maintain the number 0.7 as the constant proton fraction. 
As a significant amount of proton flux is consumed during the thermonuclear 
reaction, another situation may arise, where the proton fraction decreases grad- 
ually with time. Such a situation is given in Model-IV where the proton fraction 
decreases to 0.16 after 100 seconds with temperature-density profile same as that 
of the Model-Ill. It is evident from Fig. 4 that abundance peaks in all cases are 
obtained at mass values of 72, 76 and 80, as a result of the existence of waiting 
point nuclei 72 Kr, 76 Sr and 80 Zr, respectively. Other peaks (for example, peaks 
at mass 85, 93 and 95 for isotopes 85 Mo, 93 Pd, 95 Cd, respectively) suggest that 
the rp process flux gets accumulated at these points due to very small positive 
or negative proton separation energies of those isotopes. We have considered 
the mass number in each case for which the flux amount drops below 1% of the 
initial flux. It is evident from Fig. 4 that for Model-I, A m 93 is the region 
above which the rp process flux fall below the range of our interest, whereas for 
Model-II, the region is around A ps 95. Other abundance peaks are observed 
at masses A =97 and 101 for the isotopes 97 Cd and 101 Sn respectively, though, 
the fractions of the total flux accumulated at those isotopes are less than 0.1% 
of the initial value. In case of Models III and IV, the locations to be stud- 
ied are around A « 93 and 91 respectively. From above observations, it can 
be concluded that, for various density-temperature profiles, the rp process flux 
falls below a significant amount near mass 90-95. The observations from Fig. 4 
suggest that the end points of the rp process has a rather weak dependence on 
different X-ray burster models. 

4 Summary 

The location above which the rp process flux falls below an insignificant amount 
is calculated using the microscopic optical model utilizing the densities from the 
RMF approach and with a new phenomenological mass formula. Present result 
is compared with result obtained from another existing work[5] and the reason 
for these indifferences between the results are discussed. Our results do not 
significantly depend on the mass models. For different X-ray burster models, 
endpoints are calculated. 
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Figure 1: rp process path for 1.2 GK 




2sNi 



28 30 32 

Figure 2: rp process path for 1.5 GK 
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Figure 3: Proton capture rates (cm 3 mol 1 sec 1 ) vs temperature(GK). See text 
for details. 
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Figure 4: Relative abundances vs mass number for different density-temperature 
profiles of various X-ray burster after 100 seconds. (Insets) Model-II: Curve for 
temperature in GK vs time in seconds of the X-ray burst from Schatz et al. [2] ; 
Model-Ill and IV: Temperature in GK^fbft) and density p in gram/ cm 3 (right) 
with time in seconds from Illiadisll?]. See text for details. 



